† Corresponding author. E-mail:
Project supported by the National Natural Science Foundation of China (Grant Nos. 11834005 and 11674285).
A general approach is proposed to the quantum Rabi model and its several variants within the extended coherent states. The solutions to all these models including the anisotropy and the nonlinear Stark coupling are then obtained in an unified way. The essential characteristics such as the possible first-order phase transition can be detected analytically. This approach can be easily applied to the recent experiments with various tunable parameters without much additional effort, so it should be very helpful to the analysis of the experimental data.
The basic interaction between a two-level atom and a classical light field was described by the Rabi model many years ago.[1] Its fully quantized version in the rotating wave approximation (RWA) was introduced by Jaynes and Cummings later.[2] The Jaynes–Cummings model can be easily solved due to the conserved excitations of the atom and the photonic number. It has been widely used in the quantum optics, because the basic physics explored in the Jaynes–Cummings model alone can be realized and observed in the earlier experiments due to the extremely weak coupling between the atom and the cavity, such as Rabi oscillations, collapses and revivals of quantum state populations, quadrature squeezing, and photon anti-bunching.[3,4]
However, the situation has changed in the past decade. In many advanced solid devices, such as the circuit quantum electrodynamics (QED) system, two-dimensional electron gases, and trapped ions, the ultrastrong coupling[5,6] and even deep strong coupling[7,8] between the artificial atom and resonators have been accessed, and the RWA is demonstrated to be invalid.[5] On the other hand, the two-level system appearing in the quantum Rabi model (QRM) and its variants is a qubit, which is the building block of quantum information technologies with the ultimate goal to realize quantum algorithms and quantum computations. Just motivated by the experimental advances and potential applications in quantum information technologies, the QRM in which the RWA is not made has attracted extensive attention.[9–20] For more complete review, please refer to Refs. [21–23].
The QRM continues to inspire exciting developments in both experiments and theories recently. The anisotropic QRM[24–26] was motivated by the recent experimental progress.[27–29] It can be mapped onto the model describing a two-dimensional electron gas with Rashba (rotating wave coupling relevant) and Dresselhaus (counter rotating-wave coupling dependent) spin–orbit couplings subject to a perpendicular magnetic field.[27] These couplings can be tuned by an applied electric and magnetic field, allowing the exploration of the whole parameter space of the model. This model can directly emerge in both cavity QED[28] and circuit QED.[29] For example, in Ref. [30], a realization of the anisotropic QRM based on resonant Raman transitions in an atom interacting with a high finesse optical cavity mode was proposed. On the other hand, Grimsmo and Parkins proposed a novel scheme by adding a nonlinear coupling term to the QRM Hamiltonian.[31] This nonlinear coupling term has been discussed in the quantum optics literature under the name of dynamical Stark shift, a quantum version of the Bloch–Siegert shift, so it was later named the quantum Rabi–Stark model (RSM).[32] This model also attracted much attention in recent years.[33–35] Recently, the anisotropic Dicke model with the Stark coupling terms, which can be called as the anisotropic Dicke–Stark model, was demonstrated via cavity assisted Raman transitions in a configuration using counterpropagating laser beams.[36] For one atom case, it is just the anisotropic RSM.
With the progress on various extensions of the QRM relevant experiments, different approaches have been developed to solve various QRMs.[37] In this work, we will introduce a generical approach to solve the QRM and its several variants, such as the anisotropic QRM, RSM, and anisotropic RSM, in an unified way.
The paper is structured as follows. In Section
The general QRM can be described as follows:
![]() |
Fortunately, associated with this very general Hamiltonian is still the conserved parity
First, we employ the following transformation:
![]() |
![]() |
![](cpb_29_2_020302/cpb_29_2_020302_ieqn5.gif)
![](cpb_29_2_020302/cpb_29_2_020302_ieqn6.gif)
Next, we expand the eigenfunction as
![]() |
By the Schrödinger equation, we have the same following equation for both up and down levels due to the conserved parity:
![]() |
![]() |
![]() |
![]() |
For any real physical systems, the wavefunction should be analytic and well defined, so the higher order coefficients should be vanishingly small. Strictly speaking, cm→∞ → 0. But in the real calculation, one cannot set m to infinite as the numerical diagonalization proceeds in the truncated Hilbert space. Thus we need truncate the summation by setting cmtr+1 = 0. cmtr + 1 is actually a polynomial in α, also in the energy E if using relation (
Most importantly, for arbitrary model parameters, the saturation calculation can be arrived at if mtr is large enough, which results in exact solutions. It is very crucial to obtain the real roots of Eq. (
By this method, we can further calculate the energy spectra for the general model (
![]() | Fig. 2. Energy spectra for Δ = 0.5, U = 0.5, r = 0.5 by both present approach (solid lines) and the numerical ED (black circles). |
In the studies of the original simplest QRM, the analytical closed-form solution is very interesting and helpful.[39,40] Now, in our unified framework, we can also present the closed-form solutions to this general model by terminating the summation in Eq. (
![]() |
![]() |
![]() |
![]() |
Also considering the most complicate case, i.e., nonzero U and r ≠ 1, we plot the above two energies Eqs. (
![]() | Fig. 3. A few low levels for Δ = 0.5, U = 0.5, r = 0.1 by both the closed form (filled circles) and the exact solutions by zeros of Eq. ( |
In the recent experiments, many solid-state devices can be operated in the ultra-strong cavity–cavity coupling regime, i.e., g < 0.2. From Fig.
We have proposed a unified approach to the quantum Rabi model and its several variants. A polynomial equation with a single variable is derived by tunable extended coherent states for the most general model, which is more complicate than the often studied models. The solutions to this polynomial equation recover exactly all eigenvalues and eigenfunctions of the models for arbitrary model parameters. Closed-form solutions are also given for the first two levels, which are very accurate up to the ultrastrong coupling regime. The first-order quantum phase transitions for some parameters are observed in the unified exact solutions, and also appear in the closed-form ones. Further physics phenomena in the presence of anisotropy and nonlinear Stark coupling will be explored within this unified scheme in the near future.
Finally, we like to point out that this unified method should be very helpful to analyze the experimental data in the recent circuit QED systems from the weak coupling, ultra-strong coupling, and deep-strong coupling regimes. For some experimental device, the underlying interactions information such as anisotropy and nonlinearity are not very clear, and thus the specified model is unknown in priori. By using the present universal approach, one can easily fit the experiential data to the theoretical results based on the generalized model where all possible interactions are taken into account, and finally detect the real interactions in the experiments.
[1] | |
[2] | |
[3] | |
[4] | |
[5] | |
[6] | |
[7] | |
[8] | |
[9] | |
[10] | |
[11] | |
[12] | |
[13] | |
[14] | |
[15] | |
[16] | |
[17] | |
[18] | |
[19] | |
[20] | |
[21] | |
[22] | |
[23] | |
[24] | |
[25] | |
[26] | |
[27] | |
[28] | |
[29] | |
[30] | |
[31] | |
[32] | |
[33] | |
[34] | |
[35] | |
[36] | |
[37] | |
[38] | |
[39] | |
[40] |